    IOobject UMeanHeader
    (
        "UMean",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ
    );

    if (!UMeanHeader.typeHeaderOk<volVectorField>(true))
    {
        Info<< "    No UMean field" << endl;
        continue;
    }

    volVectorField UMean
    (
        UMeanHeader,
        mesh
    );

    volSymmTensorField UPrime2Mean
    (
        IOobject
        (
            "UPrime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );
    volScalarField Rxx(UPrime2Mean.component(symmTensor::XX));
    volScalarField Ryy(UPrime2Mean.component(symmTensor::YY));
    volScalarField Rzz(UPrime2Mean.component(symmTensor::ZZ));
    volScalarField Rxy(UPrime2Mean.component(symmTensor::XY));

    volScalarField pPrime2Mean
    (
        IOobject
        (
            "pPrime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );

// instane T
    volScalarField T_0_001
    (
        IOobject
        (
            "T_0_001",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField T_0_005
    (
        IOobject
        (
            "T_0_005",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField T_0_025
    (
        IOobject
        (
            "T_0_025",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField T_0_05
    (
        IOobject
        (
            "T_0_05",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField T_0_125
    (
        IOobject
        (
            "T_0_125",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField T_0_25
    (
        IOobject
        (
            "T_0_25",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField T_0_71
    (
        IOobject
        (
            "T_0_71",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
	volScalarField T_1
    (
        IOobject
        (
            "T_1",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
//Mean T
    volScalarField T_0_001Mean
    (
        IOobject
        (
            "T_0_001Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_005Mean
    (
        IOobject
        (
            "T_0_005Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_025Mean
    (
        IOobject
        (
            "T_0_025Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_05Mean
    (
        IOobject
        (
            "T_0_05Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_125Mean
    (
        IOobject
        (
            "T_0_125Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_25Mean
    (
        IOobject
        (
            "T_0_25Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_71Mean
    (
        IOobject
        (
            "T_0_71Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_1Mean
    (
        IOobject
        (
            "T_1Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_001Prime2Mean
    (
        IOobject
        (
            "T_0_001Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_005Prime2Mean
    (
        IOobject
        (
            "T_0_005Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_025Prime2Mean
    (
        IOobject
        (
            "T_0_025Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_05Prime2Mean
    (
        IOobject
        (
            "T_0_05Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_125Prime2Mean
    (
        IOobject
        (
            "T_0_125Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_25Prime2Mean
    (
        IOobject
        (
            "T_0_25Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_0_71Prime2Mean
    (
        IOobject
        (
            "T_0_71Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField T_1Prime2Mean
    (
        IOobject
        (
            "T_1Prime2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
   /*
    volScalarField tke_CkMean
    (
        IOobject
        (
            "tke_CkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField tke_DkMean
    (
        IOobject
        (
            "tke_DkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField tke_EpikMean
    (
        IOobject
        (
            "tke_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField tke_PkMean
    (
        IOobject
        (
            "tke_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField tke_PikMean
    (
        IOobject
        (
            "tke_PikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField tke_TkMean
    (
        IOobject
        (
            "tke_TkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    

    
    volScalarField uu_CkMean
    (
        IOobject
        (
            "uu_CkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uu_DkMean
    (
        IOobject
        (
            "uu_DkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uu_EpikMean
    (
        IOobject
        (
            "uu_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uu_PkMean
    (
        IOobject
        (
            "uu_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uu_PikMean
    (
        IOobject
        (
            "uu_PikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uu_TkMean
    (
        IOobject
        (
            "uu_TkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uu_PhikMean
    (
        IOobject
        (
            "uu_PhikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    

    
    volScalarField vv_CkMean
    (
        IOobject
        (
            "vv_CkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField vv_DkMean
    (
        IOobject
        (
            "vv_DkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField vv_EpikMean
    (
        IOobject
        (
            "vv_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField vv_PkMean
    (
        IOobject
        (
            "vv_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField vv_PikMean
    (
        IOobject
        (
            "vv_PikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField vv_TkMean
    (
        IOobject
        (
            "vv_TkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField vv_PhikMean
    (
        IOobject
        (
            "vv_PhikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );

    
 
    
    volScalarField ww_CkMean
    (
        IOobject
        (
            "ww_CkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField ww_DkMean
    (
        IOobject
        (
            "ww_DkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField ww_EpikMean
    (
        IOobject
        (
            "ww_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField ww_PkMean
    (
        IOobject
        (
            "ww_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField ww_PikMean
    (
        IOobject
        (
            "ww_PikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField ww_TkMean
    (
        IOobject
        (
            "ww_TkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField ww_PhikMean
    (
        IOobject
        (
            "ww_PhikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    

    
    volScalarField uv_CkMean
    (
        IOobject
        (
            "uv_CkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uv_DkMean
    (
        IOobject
        (
            "uv_DkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uv_EpikMean
    (
        IOobject
        (
            "uv_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uv_PkMean
    (
        IOobject
        (
            "uv_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uv_PikMean
    (
        IOobject
        (
            "uv_PikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uv_TkMean
    (
        IOobject
        (
            "uv_TkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    
    volScalarField uv_PhikMean
    (
        IOobject
        (
            "uv_PhikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );

//scalar budget
    volScalarField s_uu_Dk_1Mean
    (
        IOobject
        (
            "s_uu_Dk_1Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField s_uu_Dk_2Mean
    (
        IOobject
        (
            "s_uu_Dk_2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField s_uu_Dk_3Mean
    (
        IOobject
        (
            "s_uu_Dk_3Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField s_uu_Dk_4Mean
    (
        IOobject
        (
            "s_uu_Dk_4Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField s_uu_PkMean
    (
        IOobject
        (
            "s_uu_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );	
    volScalarField s_uu_PikMean
    (
        IOobject
        (
            "s_uu_PikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );
    volScalarField s_uu_EpikMean
    (
        IOobject
        (
            "s_uu_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );   
	
    volScalarField ss_T_1_Dk_1Mean
    (
        IOobject
        (
            "ss_T_1_Dk_1Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );   	
    volScalarField ss_T_1_Dk_2Mean
    (
        IOobject
        (
            "ss_T_1_Dk_2Mean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );   	
    volScalarField ss_T_1_EpikMean
    (
        IOobject
        (
            "ss_T_1_EpikMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );	
    volScalarField ss_T_1_PkMean
    (
        IOobject
        (
            "ss_T_1_PkMean",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT
        ),
        mesh
    );		
*/	
	
    /*
    volScalarField epsilonMean
    (
        IOobject
        (
            "epsilonMean",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );

    volScalarField nuMean
    (
        IOobject
        (
            "nuMean",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );

    volScalarField gammaDotMean
    (
        IOobject
        (
            "gammaDotMean",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );

    volScalarField nuPrime2
    (
        IOobject
        (
            "nuPrime",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );
    volScalarField nuPrime = sqrt(mag(nuPrime2 - sqr(nuMean)));


    volScalarField gammaDotPrime2
    (
        IOobject
        (
            "gammaDotPrime",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        ),
        mesh
    );
    volScalarField gammaDotPrime = sqrt(mag(gammaDotPrime2 -sqr(gammaDotMean)));
    */
